Dielectric catastrophe at the Mott transition 
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We study the Mott transition as a function of interaction strength in the half-filled Hubbard chain 
with next -nearest-neighbor hopping t' by calculating the response to an external electric field using 
the Density Matrix Renormalization Group. The electric susceptibility \ diverges when approaching 
the critical point from the insulating side. We show that the correlation length £ characterizing this 
transition is directly proportional to fluctuations of the polarization and that \ ~ £ 2 - The critical 
behavior shows that the transition is infinite-order for all t' , whether or not a spin gap is present, 
and that hyperscaling holds. 

PACS numbers: 71.10.Fd, 71.30.+h, 75.40.Cx 
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A material's response to an applied electric field char- 
acterizes whether it is a metal or an insulator. One such 
response is the static electrical conductivity at zero tem- 
perature, which is finite for a metal (or infinite for an 
ideal conductor), but vanishes for an insulator jjjj. The 
conductivity can therefore be used to probe the metal- 
insulator transition from the metallic side. A comple- 
mentary quantity is the dielectric response to an electric 
field, the electric susceptibility, \- This quantity is ex- 
pected to diverge (for a continuous transition) when the 
transition is approached from the insulating side and to 
remain infinite in the metallic phase. This phenomenon, 
termed "dielectric catastrophe" by Mott Q , has been re- 
ported for doped silicon ||. 

One possible origin of insulating behavior is the local 
Coulomb repulsion between electrons. This "Mott phe- 
nomenon" Q leads to a metal-insulator transition which 
occurs either as the electron density, n, is varied for fixed 
electron-electron interaction strength or as a function of 
interaction strength at fixed electron density |§J|. In 
this letter, we concentrate on the transition as a function 
of interaction strength for fixed electron density. Exper- 
imentally, such a transition can be induced by applying 
isostatic or chemical pressure. 

The prototype model for the Mott transition is the 
single-band Hubbard model with purely local interac- 
tion, whose Hamiltonian is 



H = - 2J Uj 



U^flifhii 



(1) 



where c it7 creates an electron of spin a at site i and 
flic, = c\ a Ci a . The hopping matrix elements ty are short- 
ranged. At half-filling, n = 1, the Hamiltonian ([!]) maps 
onto a Heisenberg model with couplings Jy = 4iy/J7 
for U — > oo and is thus insulating, while at U — 0, it 
describes a perfect metal. Therefore, a Mott transition 
must occur at some U c > 

In order to describe the dielectric response of such a 
system, one must consider the coupling to a static elec- 
tric field. Taking the field in the x-direction and neglect- 



ing overlaps between different Wannier functions (tight- 
binding approximation), we add the coupling term 



= -EX = -E^2 



XiTli 



(2) 



where X is the dipole operator (we have put q = 1), 
Xi is the ^-coordinate of the z-th site and hi measures 
the occupation of this site. Here we have assumed that 
the finite lattice has open boundary conditions, i.e. the 
connections terminate at the lattice edges. We note that 
while this is the natural definition for experiments, the 
notion of response to an applied electric field has recently 
been generalized to periodic boundary conditions JtJ . An 
applied electric field will induce a polarization at zero 
temperature given by 



P = L- d (X) 



-L~ 



dE 
dE 



(3) 



on a ci-dimensional lattice with linear dimension L, where 
the average is taken with respect to the ground state of 
the full Hamiltonian H + H ext , with corresponding en- 
ergy Eq. The zero-field susceptibility is then defined as 
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The examination of the properties of this susceptibility 
in the vicinity of the Mott metal-insulator transition is 
the principle aim of this letter. 

The susceptibility x can be related to the eigenstates 
|\& n ) of H using elementary perturbation theory, 
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where AE n is the excitation energy of the n-th eigenstate. 
(Here we have chosen the origin of the coordinate system 
so that (X) = for E = 0.) This relation immediately 
yields a useful inequality in terms of the "charge gap" , 
A (defined as the lowest excitation energy for which the 
dipole matrix element does not vanish): 
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x < 1 L- d (^o\X 2 \^ ) ■ (6) 

It is thus instructive to consider fluctuations of the po- 
larization, (\E'o|^ 2 |^'o}, which can be estimated as fol- 
lows. We expand the ground state as a series \^o) = 
J2d l^o )> where D is the number of doubly occupied 
sites ( "particles" ) . At large U the "particles" are located 
close to empty sites ("holes"). Each particle-hole pair 
represents an elementary dipole with essentially random 
orientations. Therefore our estimate is 

(* |1 2 |* ) = ^(^|1 2 |^) « (D)l 2 , (7) 

D 

where I is the average size of the dipoles. Comparing this 
result with the inequality in Eq. (^), we conclude that a 
diverging susceptibility requires either a diverging size of 
the dipoles or a vanishing charge gap or both. In one 
dimension, the quantity 

e = i(* |i 2 i*o) (8) 

is a length characterizing the insulating phase ^,|) . We 
will show below that £ is the correlation length, up to a 
dimensionless constant. 

On regular lattices, one often faces the problem that 
the Mott phenomenon, which sets in at large values of 
U due to charge blocking, is masked by the opening of 
a charge gap at much lower values of U due to antifer- 
romagnetic order induced by nesting or Umklapp pro- 
cesses. In order to control such effects, we consider here 
a model that explicitly incorporates frustration of an- 
tiferromagnetism, namely the one-dimensional Hubbard 
model with both nearest-neighbor t and next-nearest- 
neighbor t' hopping terms. We set t = 1 and consider 
only t' > here because the sign of t' is irrelevant at 
half-filling due to particle-hole symmetry. For t' = 0, 
the Bethe-Ansatz solution allows one to calculate the 
charge gap Jlfl], the charge stiffness, and the correlation 
length in the insulator |ll[] explicitly. The system is found 
to be insulating for all positive values of U. The metal- 
insulator transition occurs at U c = and is infinite order: 
the charge gap and, correspondingly, the inverse of the 
correlation length decrease exponentially as U — > + . At 
the same time, the magnetic correlations show a power- 
law decay. For t' > 0, a weak-coupling renormalization 
group analysis fl^| predicts the same behavior as long as 
there are two Fermi points: Umklapp processes lead to an 
insulating state for U > 0, while the magnetic excitation 
spectrum remains gapless. 

For t' > 0.5, there are four Fermi points in the nonin- 
teracting band structure and the picture becomes more 
complicated. In weak coupling, the lowest-order Um- 
klapp processes are marginally irrelevant and the 
system is predicted to be metallic (vanishing charge gap) 
with a spin gap. At strong coupling, the model can be 



mapped to a frustrated Heisenberg chain, which develops 
a spin gap for J'/ J ~ t' 2 > 0.2412 |b|] and incommensu- 
rate antiferromagnetic order for J' /J > 0.5. This general 
picture has been confirmed numerically [Ti|JT5[ | . For a de- 
tailed phase diagram, we refer the reader to Fig. 3 of Ref. 
p"5| . Here we will examine both the parameter regime 
with gapless magnetic excitations and U c — (f < 0.5) 
and the one with gapped spin degrees of freedom and 
U c > (f > 0.5). 

In order to numerically evaluate the electric suscep- 
tibility, Eq. (0), we use the Density Matrix Renormal- 
ization Group (DMRG) |l6|]. We apply a small electric 
field so that the system is in a linear response regime 
(typically EL = 0.001) and measure 

i 

We use the finite-size DMRG algorithm [jl6],[l7| on up to 
L = 100 sites, retaining up to 2400 states for the system 
block. This allows us to keep the sum of the discarded 
density-matrix eigenvalues to below 10 -8 . We have per- 
formed extensive tests for U — 0, a difficult case to treat 
numerically, and find that we can reproduce analytic re- 
sults to within less than one percent. The details of the 
calculations will be described more extensively elsewhere. 




FIG. 1. Electric susceptibility, \, as a function of 1/L for 
(7=1 (circles), U = 2.5 (stars), U = 4 (diamonds), U = 5.5 
(crosses), and U = 7 (squares). 

The electric susceptibility \ 1S shown in Fig. | as a 
function of the inverse system size for t' — 0.7 and a 
number of (/-values. There are two characteristically 
different behaviors: at small U , the system is metallic, 
and the susceptibility diverges with system size. A fit 
to a power law in L yields an exponent very close to 2 
(within 5%) for the small U values. For U = 0, it can be 
shown analytically that \ ~ L 2 for large L for all values 
of t'. We conjecture that such a L 2 divergence of \ is 
generic for a one-dimensional perfect metal. For larger 
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U, x extrapolates to a finite value as L — > oo. While this 
is clear for the two larger U -values in Fig. [l], care must 
be taken near the transition because the system appears 
metallic up to a length scale on the order of the corre- 
lation length which diverges at the transition. Such a 
crossover from metallic to insulating behavior is evident 
in the U — 4 curve, for which we have taken lattice sizes 
of up to L = 100 to show that x scales to a finite value, 
i.e. that the system is insulating. 

12 




FIG. 2. Electric susceptibility x°o{U,t') of the infinite-size 
system for t' = 0, 0.7, 0.8, as a function of U; the lines are 
guides to the eye. Inset: x°°{U,t' = 0.7) for U = 4.1 to 4.4 
(squares) as a function of 1/(U — U c ) , on a semilog scale. The 
line is a fit to an exponential form. 

In the insulating regime, we expect x to be analytic 
in l/L. We therefore perform finite-size scaling for large 
L using a linear fit and extrapolating to l/L = 0. The 
result, Xoo, is shown in Fig. || for t' = 0, 0.7, and 0.8 as 
a function of U. For t' — 0, the transition takes place 
at U c = 0, as discussed previously. Although we could 
not obtain a reliable finite-size extrapolation for U ^ 2 
because the correlation length becomes much larger than 
the system sizes we were able to reach, we could observe 
numerically that x ~ A~ 2 (for U < 10), where A is 
the charge gap given in Ref. pp[ . The extrapolation to 
A — > confirms that U c = 0. Alternatively, we can fit to 
the \ow-U form for A~ 2 , 



Xoo(t' = 0) 



.4 



U-U, 



exp 
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(u - u c )° 



(10) 



with the exactly known values B — At: — 12.566 . . . and 
fj = l; here the prefactor 1/(U — U c ) comes from the 
logarithmic correction. This yields U c ~ 0.058 and we 
effectively find that U c = to within error bars. A fit to 
the form without the logarithmic correction would yield 
U c ~ 0.209, which is also consistent with zero, but to 
within a larger error bar. 

It is clear from Fig. ||that the bigger t' ', the larger the 
U at which x diverges. However, one must perform care- 
ful fitting in order to accurately determine U c and the 
form of the divergence at t 1 > 0, as an analytical result 



for the charge gap exists only at t' = 0. For t' = 0.7, 
we have calculated x a t many U- -values near the transi- 
tion and have fitted to both power law, x ~ (U — U c y 
and exponential forms (Eq. |l0|), but without the loga- 
rithmic correction. The logarithmic corrections are, in 
general, non-universal, i.e. t'-dependent. Leaving these 
corrections out, as argued above, will only make the de- 
termination of U c less precise. We find that the fit to 
the power law form yields U c — 3.4, a point at which 
careful finite-size scaling of x yields a finite value of Xoc ■ 
Therefore, this U c is clearly too large. The exponential 
fit yields a ~ 1.049, B ~ 12.45 and U c ~ 2.67, a more 
reasonable value of U c . Note that the values for a and 
B are again very close to the ones obtained for t 1 = 0. 
The inset of Fig. || shows a semilog plot of Xoo versus 
1/(U — U c ) as well as the fit itself, illustrating its good 
quality. We therefore find that the exponential form, Eq. 
(|l0|), expected in an infinite-order transition, character- 
izes the transition at all t' , irrespective of whether a spin 
gap exists or whether U c is finite or zero. 

If hyperscaling is valid, there is only one relevant length 
scale £oo (the correlation length for L — > oo) in the vicin- 
ity of the quantum critical point. This length then de- 
termines the finite-size scaling of the singular part of the 
ground state energy density |is{ ] 



E s ine /L d 



C (d+2) /(£/£oo), 



(11) 



where z is the dynamic critical exponent and / a univer- 
sal scaling function. The quantity EL is an energy and 
therefore scales like £_^ Q Z . Using Eq.(^), one obtains the 
scaling behavior of the electric susceptibility 



X = L^-tC^L/^), 



(12) 



where C is a non-universal constant that depends on mi- 
croscopic details and <I> is a universal function |L9|]. The 
hyperscaling assumption also implies that $ tends to a 
finite value as L/£oo ~~ * 0. This is the region in which the 
system appears metallic and in which \ tends to scale 
like L 2 . Thus z = 1 is the only consistent value in Eq. 
(|j~2|), in agreement with exact results for t' = pi| . In 
the opposite limit, £/£°o — > oo, the system behaves as 
an insulator for all sizes and x tends to a finite value 
Xoo. The scaling form ( p^ ) with z = d = 1 thus implies 
liniz^oo $(x) — A/x 2 and 



Xoo = C A i 
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(13) 



where A is a universal constant. 

In order to confirm the scaling form Eq.jl^) for our 
model, we plot the DMRG results for x/L 2 as a func- 
tion of L/^oo in Fig. ^. The quantity ^ is obtained by 
calculating £ on finite systems using Eq. (j^) and then per- 
forming a finite-size extrapolation similar to that used to 
obtain Xoo- Notice that all L and U points for a partic- 
ular t' collapse onto the same curve, confirming hyper- 
scaling. Therefore, ^ behaves as the correlation length, 
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which we have checked by ascertaining that £oo is the 
same length (up to a constant) that characterizes the ex- 
ponential decay of the density-density correlation func- 
tion. 



%/L 2 




FIG. 3. Scaling plots of X (L, U, t') /L 2 versus L/^ x (U, t') 
in a log-log scale: t' = (left), t' = 0.8 (right). Different 
symbols correspond to different values of U. 



finite-size system depends only on the ratio £/£oo- The 
finite-size scaling of \ can then be related to a universal 
scaling function and a dynamic exponent z = 1. The 
transition is found to be infinite order and to show the 
same critical behavior whether there is a spin gap or not. 
As to the origin of this "dielectric catastrophe" , we con- 
clude, on the basis of both the inequality x < 2£/ A and 
the observed scaling Xoo ~ £m> that it involves both a 
diverging correlation length £ (linked to the unbinding of 
dipoles) and a vanishing of the charge gap A. 
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The results of the l/L extrapolation for x an d £ 
are shown in Fig. [| for three different values of t'. A 
power-law fit to Xoc = C'£2o yields j(t' = Q) ~ 1.97, 
7(i' = 0.7) ~ 2.01, 7(t' = 0.8) ~ 1.96, and C'(t' = 
0.7)/C"(f = 0) ~ 1, C'(t' = 0.8)/C'(t' = 0) ~ 0.7. This 
confirms the scaling behavior ( |l3|) . It also shows that the 
constant C depends weakly on t' . 
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FIG. 4. Electric susceptibility x°°> versus correlation 
length for different values of t' . The lines are power-law 
fits. 

In summary, our calculations for the U — t — t' chain 
at half-filling confirm that the electric susceptibility x 
(and therefore also the dielectric constant e = 1 + 47rx) 
diverge when approaching the Mott transition from the 
insulating side. The polarization fluctuations, which also 
diverge for U — > U c from above, have been found to be 
directly proportional to the correlation length £ of the 
Mott insulating phase. In agreement with the hyperscal- 
ing hypothesis, the metallic or insulating behavior of the 
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